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NATIONAL ADVISORY CCMIITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 133^ 


THE EFFECT OF HIGH VISCOSITY ON THE FLOW AROUND 
A CYLINDER AND AROUND A SPHERE*^ 

By F. Homann 


For the determination of the flow velocity one is accustomed to 
measure the impact pressure, i.e., the pressure intensity in front of 

an obstacle. In incompressible fluids the Inpact pressure is 

{, 7 } is specific weight and v, m/sec, is velocity) if the 

Influence of viscosity can be neglected. Such an Influence is appreci- 
able, however, when the Reynolds number corresponding to in 5 )act tube 
radius is under about 100, and must conseq.uently be considered, if the 
velocity determination is not to be faulty. The first Investigations 

p 

of this Influence are included in the work of Miss M. Barker . In the 
following pages, e^fperiments will be reported which determine the inten- 
sity of Impact pressure on cylinders and spheres; furthermore a theory 
of the phenomenon will be developed \dxich is in good agreement with the 
measurements . 

The research apparatus consists of an oil circulation in which the 
velocity of the oil can be varied from 0.5 centimeter per second to 
30 centimeters per second with the help of a vane-type punp lying 
entirely in the oil. A Russian bearing oil and a mixture of this with 
fuel oil is used for the measurements. Figure 1 illustrates the test 
setup. In this is Indicated: P, the puiip; a, turning vanes; 

G, strai^tener; and V, the actual test section which possesses a 
breadth of 0.148 meter, a depth of 0.15 meter calculated from the oil 
surface, and a length of 0.74 meter. It was provided with wall ports A 
in three different places. E is an entrance section for the pump; 

D, a diffuser; the immersion heater T and the cooling coil K provide 


*"Der Elnfluss grosser Zahigkeit bel der Stromung urn den Zylinder 
und urn die Kugel", ZAMM, vol. l6, no. 3^ June 1936, p- 153-164. 

^The suggestion of the present work, which was prepared in the 
Kaiser Wilhelm Institute in Gottingen, I obtained from Herr Professor 
Dr. Prandtl, to whom in this place I express my most heartfelt thanks 
for the energetic furthering of the work and the valuable suggestions 
given me for its completion. Another work in the same field is published 
in "Forshung a.d. Gebiet des Ingenieurwesens", 1936, vol. 7, no. 1. 

%I. Barker - Proc. of Royal Society, 1922, Vol. A. -101, p. 435. 
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temperature regulation. For Impact pressure measurement a cylinder 
which was provided with a port was built rigidly into the test section. 
The diameters of the cylinders used were 1 centimeter, 1-377 centimeters, 
1-953 centimeters, and 2 centimeters. The holes had a diameter of 
0.1 centimeter and 0.2 centimeter. Two corrections - one because of 
wall effect, the other because of finite size of hole (which originated 
with A. Thom and first had to be checked for the measuring range under 
consideration here) - were applied to the measurements, which are illus- 
trated in figure 7- The solid curve represents the theory. 

More precise Information on the test setup and the measurement 
technique are found in the work cited in footnote 1. 

In the case of the measurement of static pressure on a sphere, a 
sphere provided with a hole was affixed to a pitot tube, the sphere 
having one time a diameter of 0.8 centimeter, the other time 1.6 centi- 
meters. The execution of the measurements was in the same manner as in 
the case of measurement of the Barker effect. The result is shown in 
figure 8. The solid cu 2 ?ve corresponds to the theory. 

In order to arrive at a clearer picture of the viscous flow around 
a cylinder or a sphere, the case of viscous flow against a plate was 
next calculated. The differential equation appearing in the two- 

dimensional case has already been solved by Hlemenz^ and will be sketched 
once more for the sake of a better understanding of the final form. The 
solution will than be used in the flow around the cylinder. 

After this, the three-dimensional flow against a plate in a fluid 
jet will be treated, to be used on the sphere. 


VISCOUS FLOW IN THE VICINITY OF A STAGNATION POINT 
(TWO-DIMENSIONAL CASE) 


Potential- flow theory gives for the velocity components in the 
neighborhood of a stagnation point, for the case of flow perpendicular 
to a plane wall (fig. 2): 


U = -ax V = ay 


%iemenz - Dissertation, Dlngl. Polyt. Journal, V .326 (I9II), 
No. 21-26. 
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The pressure is found from the Bernoulli equation to be 
P„ - p = + v") - . 

With consideration of viscosity, Hiemenz^ makes the foimulation; 


u - -f(x) 

V = yf’(x) 

(1) 

Po-P-2|^( 

^(x) " 

(2) 


The continuity equation is fulfilled; the boundary conditions read: 

for X = 0 (that is, at the wall): u=v=0, f=f*=0 

for X = oo: V = V, f’ = a 


In equation (2) if signifies the pressure at the stagnation point, 

then From the equations of motion 


tc- 

C) C 
c-'-tt- 

<f t 


^u Su 

U + V 

_ 1 Sp' j 


a^u' 
+ ■ - 

5x 


“ “p ^ ' 

vax^ 

3y^ 

Sv Sv 

U V V 

= -i^+ v| 


S^v' 
+ 

Sx 

Sy 

p dy ' 

[dx^ 

3y^, 


one obtains as determining equations for f, v and F/ \ : 

(x) (x) 


ff» = — F* - Vf* * 
2 


m2 


- ff*« = 


Vf 


t t t 


( 3 ) 

(^) 


with the boundary conditions above. 
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Equation ( 4 ) has already been integrated by Heimenz. In order to 
make the coefficients equal to unity, he sets 

f(x) = A 0 (|) I = ax 

From comparison of the coefficients, it follows that 

A = ^ ^ 

With this, equation ( 4 ) becomes 

»0 = 1 + » » 

The new boundary conditions read 

for 1=0: 0 = 0* = 0 

for I = 00: 0’ = 1 

The behavior of 0 and the first two derivatives is shown in figure 4 . 

We need the pressure difference between the stagnation point and 
the pressure for x = For x = 00: 

0‘=1 f‘=a 0=1- 0.647 

From integration of equation ( 3 ), F is determined to be 

^ F = Vf + I f 2 (6) 

If one now forms (p^^ - p^ minus ^p^’ - p')^ as given by the Bernoulli 
equation, one obtains: 



I 


( 5 ) 
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For the stagnation streamline, for which y = 0: 

K - f) - (Po' - P')= ^ 

If one puts in for F the previously obtained value, there results 
(Pq - p) - (p,,‘ - P’) = |(2vf J + f^2 _ 

For X = 00 , f* = a; therefore, one obtains as a final formula: 

(Pq - p) - (Pq' - P') = Pva (6a) 


VISCOUS FLOW AT A STAGNATION POINT 
(ROTATIONALLY SYMMETRIC CASE) 


For the solution of the differential equation arising, all the 
expressions, such as the equations of motion, the velocity components, 
etc,, were reduced to cylindrical coordinates. 


If z, r and p are the coordinates (fig. 3)^ then corresponding 
to the two-dimensional case, there will, apply: 

(7) 


= -^(z) 


V. = I f* 


r p " (z) 


Po - p 




( 8 ) 


The continuity equation is again fulfilled; Vp = 0, since we are dealing 

with a rotationally symmetric process. These expressions stem from the 
frictionless problem of a fluid Jet against a plate, where 


and 


Vj. = ar V 2 = -2az 

. p . E|!(h. 2 , rU 
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The quantity 2az 
the viscous case. 


In the frlctlonless case Is replaced by f 


(z) 


in 


In the case at hand the equations of motion read: 










Sz 


P Sr 



1 

^ 5r 



1 ^ + 
P Sz 



S'- Sz2 / 


J 


(9) 


Substituting equations (j) and (8) In equation ( 9 ) gives: 

( 10 ) 
( 11 ) 

The boundary conditions read 

for z = 0: f = f« = 0 

for z = «>: f’ = 2a 


1 f«2 _ fpii 2 a^ + vf ' ” 

2 




-t 




fffQ = 


a*- 


F‘ - Vf* 


0 


If one finds f from equation (lO), one can there-vri-th determine F 
from equation (ll). One next substitutes the transformation 

f(z) = A0(O I = az (12) 

Into equation (10), In order to make the t-wo coefficients equal to unity. 
This yields: 

i a^A^0*^ - a%.^00* ’ = 2a^ + va^A0’*’ 


From equating the coefficients: 

1 2 2 2 1 
- a A = 2a^ = Va^A 
2 
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A = 2 ^v 



(13) 


From equation (lO) with equations ( 12 ) and equations (I3) there results 
the final differential equation: 


j3»’» +\200** 


0 *^ + 1=0 


with the boundary conditions: 


( 14 ) 




for 1=0: 0=0* =0 

for I = c»: 0’ = 1 


The differential equation (l 4 ), just as Hiemenz*, is no longer 
elementarily integrable. Its solution was obtained, accordingly, 
through a power series development from zero: 


0 = a^ + aj 


+ Sir 




n 


(15) 


By the method of undetermined coefficients, a^ can be determined'. 

Since, however, one boundary condition lies at infinity, one coefficient 
remains undetermined; and in fact it turns out to be a.2. From the 

recursion formulas 


2 3 

0 = a^^ + a^^ + a^l + + . . . 

0’ = a^ + 2a2l + + . . . 

0** = 2ag + 2 X 3^2 1 + 3 X 4aj^|^ + . . . 


result as coefficients: 
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ao = = 0 

a-2 = for the present undetermined 

a3 = -0.166667 
= 0 
3.^ = 0 

= 0.555556 X 10"^ a.2 

arj = -0.396825 X 10"3 

ag = 0 

a^ = -0.440917 X 10“3 0,^2 

^10 = 0. 793651 X 10-^ ag 

^11 = -0.360750 X 10-5 

^12 = 0.374111 X 10-^ ag3 
aj^g = -0.114597 X 10-^ ag2 
a^i^ = 0.115735 X 10-5 ag 

ai5 = -0.301482 X 10"5 ag^ - 0.385784 X 10" 

= 0.134896 X 10"5 ag3 
= -0.211005 X 10-^ ag2 

a^8 = 0.224141 x 10"^ ag5 + 0. 157758 x 10"^ ag 
^15 = -0.135546 X 10-^ ag^ - 0.415153 X 10"9 
agQ = O.316633 X 10"7 ag3 

agi = -0.152798 X lO-T ag^ - 0.295658 X 10“^ ag2 
^22 = 0.119505 X 10-7 ag5 + 0.199390 X 10"9 ag 

^23 " -0.371665 X 10~® ag^ - 0.433457 X 10-11 
^24 " 0.956242 X 10"5 ag7 + 0.554360 X 10-9 ag3 

ag5 = -0.943031 X 10-9 ag6 - 0.462914 x 10"11 ag2 


J 
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In order now to be able to determine 3 - 2 ^, a second series develop- 
ment from Infinity was set up, which was adjusted to the boundary condi- 
tion for 0 at infinity. To this end one sets 

0 = 0Q ^1 

in which 0^^ corresponds to a small quantity, which one can neglect -in 
the following expressions when it appears squared. 0^ is the solution 
for I = 00. 


0’ = 0o* + 01* = 1 + 01* 


since for | 


0Q» = 0‘ = 1 0»« = 0^«« 0«»« = 02* ** 


The boundary condition reads 


for 5 = 00: 0 « z= 0 


(16a) 


Furthermore, 0 q = 


The integration constant is omitted, since in the following calculation 
it comes in again automatically. 

If one substitutes the above values into equation (lO), one obtains 

0i'" + 2(^o0i” + 0A”) - (0o'^ + A'0i' + 1 = 0 


Or if one neglects the squared terms in 0 ^: 

03 _“' + 2 | 0 ^»* 


- 20 ^* = 0 


( 18 ) 
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To solve this differential equation one sets 


$ = 03 ^* 


With this, equation (I8) gives 


$1 1 4. 2|$i _ 20 = 0 


(19) 


A special, not identically vanishing solution of equation (19) is: 


= I 


If $2 is an additional solution, then 


- 


25 dl 


® 


- $2 = 


O2’ - i O2 - i e“^ = 0 


This equation is directly solvable. Its solution is 


= 1/^ e-^^ dq 

O 00 


The general solution of equation (I9) is then: 

2 

e ^ dq 
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Since 


for I = 00 , ^^’ = <& = 0, then = 0. Therefore 


= $ = C2S J' dri = C2I - -i e“^^ - 2 e“^^ 


di] 




and since (H-, = j dl 

CO 


p I 2 pn 2 

02^ = + C 2 / -e”^ - 2 t] / e”*^ du 


dTi 


= C 3 - C 2 


2 r ^ 2 

e“^ dri - 2C2 / r\6.-r\ e“^ 

J CD J CD 


du 


00 ^ 00 


The double integral becomes, according to Blasius^, 


2C2 J Ti di^ e-*^^ du = 2C2 I ^ ^ J di] + 


With this, equation (20) becomes: 


01 - C 


3 _ 1 


le-t". ( 1 ^ 1")/^ eV 


dq 


Now: 


ni 2 „2 P“ „2 „2 

2 / e"^ dT| = 2 / e"’^ dq - 2 / e“^ drj = 2 / e“^ dq 

'-'n '-'r» 


- ^/iF 


Vif^ 




^Blasius - Dissertation, Z.F.M.u. Phys. I 90 S, p. 1, 


( 20 ) 



( 21 ) 


>( 22 ) 
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If one substitutes equation (22) In equation (21), one can calculate ^ 
polntwise, since 





e-T 


dq 


is tabulated. Therefore 


01 - C3 




e“'n dtj - 1 1 




1 ^ 


0i” 

C2 






Herewith 0, 0’, and 0*’ are determined for the development at 

infinity, as a comparison with equations (16) and equation (l 6 a) shows. 

In both developments a 2 , ^3 appear as unknowns. If 

one now combines both solutions at the point | = |q, and determines 

that the value of the function and the first two derivatives of the 
series development at zero are equal to the corresponding values that 
one obtains from the development at infinity, then three determining 



equations for C2 and result, from which the unknowns can be detennined. Therefore: 



Iq was chosen 1 . 8 , since 0 ’’, from which 0 is built up, can then be determined accurately 
to 0 . 002 , on account of the alternating signs of the power series development. 

The solution of the determining equations gives 

a^ = 0.658619, C2 = 2.16492, = -0.557611 


H 

U) 
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Thereby a .2 is determined accurately to at least five places. For the 
coefficients of the power series this yields: 


a2 = 0.658619 

^16 = 

a^ = -0. 166667 

II 

a^ = 0.365900 X 10 "^ 

00 

II 

a.^ = -0.396825 X 10 " 3 

^19 “ 

a^ = -0.191261 X 10 “ 3 

II 

0 

OJ 

ni 

a^Q = 0.522714 X 10 “^ 

^21 

a^^ = -0.360750 X 10“5 

»22 = 

a^2 = 0.106882 X 10"^ 

^23 " 

= -0.497098 X 10"5 

ag 4 = 

a^4 = 0.762253 X 10“ 6 

"25 = 

a^^ = -0.605859 X 10"6 



0.385391 X 10-6 

-0.958673 X 10 “ 7 
0.381678 X 10“'^ 

-0.259200 X 10 “ 
0.904605 X 10 "® 

-0.252966 X 10"® 
0.161233 X 10“6 
-0.703675 X 10"9 
0.209783 X 10"^ 
-0.970520 X 10" ^6 


a 

o 


= a. = a,. = 



^8 " 


0 


The values for , and 0* ’ are shown more accurately how- 

ever, in Table I. In this case, ’ is calculated accurately to two 
decimal places, 0* to two, and 0 to three. With this the differ- 
ential equation (l4) is solved. 

From integration of equation (ll), one obtains 

^ F = Vf* + i f^ = 2av(0‘ + <0) 

■ < 


(23) 
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As in the plane case, one uses again the pressure difference 
between the stagnation pressure and the pressure for z = ». For 1=0, 
equation ( 23 ) is equal to zero; for | = oo 

0* = 1 f» = 2a 0=5- 0.557611 

If one now forms again ^p^ - pj minus (Pq' “ P*)^ given by 
the Bernoulli equation, one obtains: 

(Pq - P) - (Pq’ - P') = p(vfoo' + I foo^ - I foo^) = PVf„' 

As a final formula one obtains 

(Pq - p) - (Pq' - P*) = 2pva (24) 

STAGNATION PRESSURE ON A CYLINDER 


For the stagnation streamline the Navier-Stokes differential equa- 
tion gives 


+ (25) 

6 x P ax 3^2 ^y2 


Figure 12 shows the variation of u on this streamline. The different 
behavior of u at the stagnation point from potential flow is explained 
by the influence of viscosity. If one integrates between the boundary R 
and 00 , one gets 




V 



ax 



~ dx 


Figure I 3 shows that 
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likewlse Uj^ = 0, and we want to identify pj^ with of the pre- 

ceding calculation. One obtains, therefore 

u 2 

(Pq - p) - = PV / ^ ix ( 26 ) 

^00 Sy 



^ dx 


is calculated approximately in that for u the value corresponding to 
the potential flow is put in. The contribution of the boundary layer 
to the integral is, in the case of not too small Reynolds number, small 
in comparison. As potential function $ of the flow around the cylinder, 
one obtains 



and with it: 

S^u ^ u / ^ ^ 8x%^ 8y^R^ 48x^y%^ \ 

^ " r6 r6 ■ ,8 / 

For y = 0, therefore, along the stagnation streamline: 




PV 


.,R 2 

2-adx 

» 5y^ 


2pvUc 

R 


Herewith equation (26) gives 


(Po - P) = 


pUq^ 2pvU^ 


2 


+ 


R 


(27a) 
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or 


pUo2/2 


( 27 ) 


If one substitutes 7 = pg, then formula (27) reads 

= 1 + ( 27 b) 

Re 


where g = 9«8l meters per second^. 

In order to be able to accomplish a comparison of test results with 
theoiy, the "displacement thickness" (see Toll mi en; Hdb. d Experimental- 
physikj V. 4, 1st part, p. 262 , "Grenzschichttheorie (Boundary Layer 
Theory)") on the cylinder must yet be considered in the calculation. 
Solution of the differential equation ( 5 ) yields (fig. 5) 5 


Pq - P 

7UoV2g 


I* = (XX* = o8* = 0.6^7 


where 8* is the displacement thickness. Therefore 


0.647 = 


( 28 ) 


If one compares the flow in the region nearest the stagnation point for 
the cylinder and for the flow against a plate, one obtains from equa- 
tions ( 27 a) and (6a) 


pva 


2pVUo 

R ^ 



If one substitutes this value in equation (28), one obtains for the 
displacement thickness 


0.647 



5 * = 


8 ^ = 0.457 


The dependence of 8*/R 


on 


Re is shown in figure 6. 
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Nov since in the test results Re is formed from the cylinder 
radius R, the actual effective radius is therefore (R + 5 *), and equa- 
tion (27b) is altered to: 


Pq - P 

7UoV2g 


1 + 


Uo(R + 5*) 


With this one obtains as a final rule for the stagnation pressure on the 
cylinder 


Pn - P 

-P - = 1 + 




Re + 0.457i/Se 


(29) 


In figure 7 the solid curve again gives the theory, ^ich agrees very 
well with the practice. 


STAGNATION PRESSURE ON A SPHERE 


Corresponding to a cylinder, for the stagnation streamline of a 
sphere 


8 u 1 8 p 

u — + — — = 

8x P dx 


V 


d^u /d^ d^\ 

+ v( + 

dx^ \dy^ dz^/ 


If one integrates again over x 


from 00 to R: 




+ pv 




since 



dx^ 


dx = 0 , 


U 


R 


0 


(30) 


The integral on the ri^t side of equation (30) one again solves by 
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substituting for u the value for the potential flow. The potential 
function of this flow Is 




For potential flow It Is further true that 

B^u S^u S^u _ ^ 

Therefore 



R 


If one substitutes In the last formula the value for du/^x given by j), 
one obtains 



Herewith equation (30) becomes 


P - P 
o 


pUr 


3pVUc 

R 


or 



or with 7 = pg 


Pq - P 


1 


+ 


6 

Re 


(31a) 


(31) 


(31b) 
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If here one also puts the displacement thickness into the calculation, 
one gets (since in the rotationally symmetric case = 0.5576): 


0.5576 = 6*1^ 


Comparison of equation (24) and equation (jla) yields 


Eo 

2 R 


0.5576 = 


^ \f^ 



6 * 


(32) 




It appears that the displacement thickness I'or a sphere and a cylinder 

are equal within i percent, although the displacement thickness in the 

case of plane flow against a plate is different from the corresponding 
three-dimensional flow. 

If one considers the displacement thickness in equation (31b), one 
obtains as a final stagnation pressure formula for a sphere 


£o_UL.i, 6 

/Uo^/2g Re + 0.455^^ 


(33) 


The solid curve in figure 8 corresponds to the theory; the agreement 
with test results is again satisfactory. 

From the final stagnation pressure formula the dependence of the 
numerical factor c on Re can be determined, il’ one sets 


Pq - P 
7UoV2e 


For the sphere there results 


6 Re 

Re + 0.455\^ 


(3U) 
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From Stokes' calculation one obtains for small Re: € = 3 . In 

figure 9 is drawn log Re as abscissa, e as ordinate. In the region 
from about Re = 0.1 to Re = 1 the course of e Is essentially dif- 
ferent, since Stokes' law describes an approximation for very small 
Reynolds number and the above law is an approximation for large Reynolds 
number. 

For the cylinder one obtains In the same fashion 


4Re 

Re + 0 . 457 ^^ 


(35) 


According to Lamb^, for small Re, for which the validity of the formula 
extends to about Re = 0.5: 

6 = Re ( 36 ) 

1.309 - ZnRe 


In figure 10 is again shown the dependence of g on log Re. Within 
the accuracy of measurement the test results here also confirm equa- 
tion ( 29 ). 

With the help of the flow against a plate it is now also possible 
to establish approximately the course of u, Bu/Sx, and from this p, 
on the stagnation streamline. A single curve was assumed in which, 
inside the displacement thickness &*, the magnitudes as given by the 
flow against a plate were used. From the displacement thickness on, 
which had a value of 0.0455 in the foregoing case for Re = 100, the 
potential flow was calculated. To explain the transition from viscous 
to potential flow, I would like to go through the calculation of u as 
an example. The solution of the viscous problem has as asymptote 

the tangent to the curve U 2 /u^, which was determined from potential 

theory, at the point &* = 0.0455 centimeter. In figure 11 this tangent 
is labelled t. The difference k between the asymptote t and Uj^ 

at the point x^ gives the deviation of viscous flow from potential 
flow at this point. Therefore to the value u^ at the point x^ was 

added the proper k. With the help of this procedure one obtains point- 
wise the transition from uj_ to U 2 . 


^amb - "Hydrodynamics" (2nd Edition 1931; German Edition by 
E. Kelly, p. 696 , par. 343)- 
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Su/Sx was determined correspondingly; the pressure p was found 
from the equations of motion to be. In the case of the sphere: 


Pq - P = 6 

7Uo^/2g ^00 200 Sx Uo(R + x)^ 


( 37 ) 


Instead of 6 In the last term of the preceding equation (37)^ in 
the case of the cylinder one gets the factor 4. Figures 12, 13, and l4 
are the results; by way of comparison the corresponding curves for the 
cylinder and the sphere are shown on one sheet. The curves are true, 

as already said, for Re = 100, In which R = 0.01 meter; Uq = 1 meter; 

o 

V = 0.0001 kilogram X second per meter was assumed. 

Naturally the last curves give only an approximation, which can be 
made essentially better throu^ a second approximation; yet this task 
In the framework of the foregoing work would lead too far. 


SUMMARY 


In the foregoing work the stagnation pressure Increase on cylinders 
and spheres brou^t about through the Influence of large viscosity, was 
reported on. 

For the three-dimensional problem, hence the flow around a sphere, 
a differential equation was set up which corresponded to that of Helmenz, 
who had already solved the two-dimensional case. The solution was 
ascertained likewise through an approximate method. The solutions for 
the two- and for the three-dimensional case were used for the flow 
around the cylinder and sphere respectively; the formulas so obtained 
for the stagnation pressure Increase stood in good agreement with the 
reported test results. Finally, a procedure to determine the velocity 
and pressure variation, as well as the variation of Su/dx on the 
stagnation streamline was shown and used on the practical case of 
Re = 100. 


Translated by D. C. Ipsen 
University of California 
Berkeley, California 
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TABLE I 


1 

0 

0’ 

r 

1 

0 



0 

0 

0 

1.3172 

1.4 

0.8546 

0.9476 

0.1697 

.1 

. 00614 - 

.1267 

1.2172 

1.5 

.9502 

.9635 

.1301 

.2 

.0250 

.2434 

1.1173 

1.6 

1.0472 

.9762 

.0895 

.3 

.051^8 

.3502 

1.0181 

1.7 

1.1453 

.9863 

.0622 

.k 

.097I4 

.4471 

.9200 

1.8 

1.2424 

.9905 

. 04 l 8 

.5 

.1439 

.5343 

.8235 

1.9 

1.3436 

.9935 

.0276 

.6 

.2012 

.6129 

.7298 

2.0 

1.4430 

.9962 

.0180 

.7 
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Figure 4.- 


0, 0', 0". The curves drawn out illustrate the two-dimensional 
solution, those not drawn out the three-dimensional. 





Figure 5, 



Figure 6, - Momentum thickness on a cylinder 
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Figure 9.- e for sphere. 
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Figure 10.- e for cylinder. 
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Figure 11.- 


u 



Illustration of interpolation for transition from viscous to 
potential flow. 
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Figure 12.- Velocity variation on stagnation streamline; cylinder: Curve I; 

sphere: Curve n. Re = 100. 
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Figure 13.- Variation of on stagnation streamline; cylinder: Curve I; 

sphere: Curve It. Re = 100. 



Figure 14.- Pressure variation on stagnation streamline; cylinder: Curve I; 

sphere: Curve n. Re =100. 
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